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ABSTRACT 

Models of the solar transition region made from lines other than those of helium cannot 
account for the strength of the helium lines. However, the collisional excitation rates 
of the helium resonance lines are unusually sensitive to the energy of the exciting elec- 
trons. Non-thermal motions in the transition region could drive slowly-ionizing helium 
ions rapidly through the steep temperature gradient, exposing them to excitation by 
electrons characteristic of higher temperatures than those describing their ionization 
state. We present the results of calculations which use a more physical representation 
of the lifetimes of the ground states of He I and He II than was adopted in earlier work 
on this process. New emission measure distributions are used to calculate the tem- 
perature variation with height. The results show that non-thermal motions can lead 
to enhancements of the He I and He II resonance line intensities by factors that are 
comparable with those required. Excitation by non-Maxwellian electron distributions 
would reduce the effects of non-thermal transport. The effects of non-thermal motions 
are more consistent with the observed spatial distribution of helium emission than 
are those of excitation by non-Maxwellian electron distributions alone. In particular, 
they account better for the observed line intensity ratio 7(537.0 A)//(584.3 A), and 
its variation with location. 

Key words: line: formation - Sun: atmospheric motions - Sun: transition region - 
Sun: UV radiation. 



1 INTRODUCTION 

The resonance lines of He I (at 584.3 A) and He II (at 
303.8 A) show unusual behaviour when compared with other 
strong emission lines in solar euv spectra. Reviews of ob- 
servations of the helium resonance lines and at temp ts to 
model their formation can b e found in Hammer (1997) and 
Macpherson & Jordan ( |l999| ) (hereafter MJ99), but the de- 
tails relevant to the present work are summarized here. 

Early observations (e.g. Tousey 1967) showed that the 
helium resonance lines are reduced significantly in intensity 
in coronal holes compared with the average quiet Sun, while 
other lines formed at similar temperatures show only very 
small reductions in intensity. In the quiet Sun, the helium 
resonance line intensities are too large to be reproduced 
by emission measure distributions that account for other 
transition region (TR) lines, a problem identified by Jordan 
(1975), who found disagreements by factors of 15 for He I 
and 6 for He II. Corresponding factors of at least 10 and 13, 
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respectively, were found by MJ99. Radiative transfer calcu- 
lations have also been unable to reproduce the observed res- 
onance line intensities without invoking either a 'plateau' in 
the model atmosphere (see e.g. Vernazza, Avrett & Loeser 
1981; Andretta & Jones 1997), which leads to other lines 



formed in the lower TR being over-estimated (smith 2000), 
or some departure from equilibrium (e.g. Fontenla, Avrett 
& Loeser 1993; Avrett 1999). These results have been taken 
to imply that some process preferentially enhances the he- 
lium resonance line intensities with respect to other lines 
formed in the TR in the quiet Sun, and that the enhance- 
ment is reduced in coronal holes. Recent results from the 
Solar and Heliospheric Observatory (SOHO) show that the 
helium resonance line intensities are reduced by factors of 
1.5-2.0 (Peter 1999; Jordan, Macpherson & Smith 2001) in 
coronal holes. 

Many possible enhancement processes have been inves- 
tigated, but a completely convincing explanation has not 
yet been found. The effects of photoionization by coro- 
nal radiation have been suggested as a na tural expla nation 
of the coronal hole/quiet Sun contrast (Zirin 1975), and 
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photoionization-recombination (PR) appears to be impor- 
tant in the formation of some lines in the helium spectrum 
(He II 1640.4 A - Wahlstr0m & Carlsson 1994, He I 10830 
A - Andretta & Jones 1997). Although PR probably con- 
tributes to the formation of the He I 584. 3-A line, evidence 



(Milkey 1975; Andretta & Jones 1997; Andretta et al. 199£) 



suggests that the He II 303. 8- A line and, to a lesser extent 
the 584.3-A line, are formed principally by collisional exci- 
tation. If this is the case, then because the helium resonance 
lines have unusually large values of W/kTe, where W is the 
excitation energy, their collisional contribution functions are 
sensitive to excitation by suprathermal electrons. Any pro- 
cess exposing helium ions to larger populations of suprather- 
mal electrons than in equilibrium will tend to increase the 
collisional excitation rates of the helium lines, while lines 
with smaller W/kTe will be relatively unaffected. 

Two processes by which this might occur were suggested 
by Jordan (1975,1980): the transport of high energy elec- 
trons from the upper TR or corona to the lower TR, or the 
transport of He atoms and ions by turbulent motions to re- 
gions of higher electron temperature. Both processes would 
be expected to depend on the magnitude of the tempera- 
ture gradient, which could explain the coro nal hole/quiet 
Sun contrast, since Munro & Withbroe ( 1972[ ) found dT/dh 
to be an order of magnitude smaller inside coronal holes. 

Shoub (1983) investigated the former process in his 



study of the shape of the electron velocity distribution func- 
tion (EVDF) in the transition region. His calculations sug- 
gested that EVDFs in the lower TR should have more heav- 
ily populated suprathermal tails than the local Maxwellian 
distribution. He found that this could lead to enhanced col- 
lisional excitation and ionization in helium, but did not 
calculate intensities to compare with observations. Ander- 
son, Raymond & Ballegooijen (1996) did calculate intensi- 
ties for the He II resonance line under excitation by several 
different non-Maxwellian EVDFs, and found enhancements 
over Maxwellian collisional excitation, but their calculations 
could not reproduce the observed intensities of all the lines 
they studied simultaneously. 

Radiative transfer calculations of the helium resonance 
line intensities using EVDFs approximating those derived by 
Shoub (1982,1983) have been performed, and are reported in 
a companion paper (Smith, in preparation - hereafter paper 
II). The calculations suggest that the He I and He II reso- 
nance line intensities could be increased by non-Maxwellian 
collisional excitation, but that this process would produce 
signatures in the line ratios that contradict observations. 

The second process is the subject of this paper. The 
temperature gradients in the solar TR are such that veloci- 
ties of order 10 km s _1 could carry material into regions of 
significantly higher electron temperature. Such mixing could 
be caused by non-thermal motions in the TR that can be 
inferred from observations of the excess widths of optically 
thin lines (e.g. Berger, Bruner & Stevens 1970; Boland et al. 
1973; Doschek et al. 1976; Chae, Schuhle & Lemaire 1998. 



See also Section 2.3). Helium ions would reach statistical 



equilibrium at the new local temperature much more slowly 
than the bulk of the material owing to their long excitation 
and ionization times. The helium resonance lines would be 
excited collisionally at higher temperatures and hence with 
greater rates than would be the case in statistical equilib- 
rium at the temperatures determining the ground state pop- 



ulations, while ot her tr ansition region lines are relatively un- 
affected. Jordan ( |1980| ) investigated the effects of the process 
on the He II 303. 8-A line, and a more exten sive study has 
been made recently by Andretta et al. (200C), who termed 
the process 'velocity redistribution.' 

This work represents the application of similar methods 
to more specific examples of temperature gradients, which 
are no longer assumed to be constant over the path of the 
moving clump of plasma. A new treatment is used for the 
mean lifetime of the helium ground state, and the analysis is 
extended to the cases of the He I 584.3-A and 537.0-A lines 
in a first approximation. 

In Section |^ earlier work on non-thermal transport of 
helium is reviewed, and the methods and atmospheric pa- 
rameters used here to calculate the enhancements of the 
helium line intensities are introduced. In Section ^ the cal- 
culations of the velocity redistribution enhancement factors 
for the Hen 303.8-A line and the He I 584.3-A and 537.0-A 
lines are described in detail. The possible effects of non- 
Maxwellian electron distributions in the transition region 
on the calculations are assessed. In Section [| the results 
of the calculations are discussed. The derived enhancement 
factors are compared with the results of MJ99, to determine 
whether the process can explain quantitatively the anoma- 
lously large helium line intensities. More qualitative compar- 
isons are made with observations of the spatial variations of 
the helium line intensities in the quiet Sun with respect to 
each other and to other transition region lines. Comparisons 
are also made with the results of a separate investigation into 
the possibility that the helium resonance line intensities are 
enhanced by excitation by non-local suprathermal electrons 
(Paper II). 



2 METHODS OF CALCULATION 
2.1 Earlier work 

Jordan (1980) investigated the enhancement that would oc- 
cur in the He II 303.8-A line intensity if the line were col- 
lisionally excited at a higher temperature than that appro- 
priate to conditions of statistical equilibrium. Starting at an 
initial temperature of Tj = 8 x 10 4 K, she assumed that 
He II ions are carried up the temperature gradient at the lo- 
cal non-thermal (macroturbulent) velocity derived from the 
excess widths observed in lines formed near the initial tem- 
perature. She wrote the excitation time of the ion's ground 
state as the reciprocal of the collisional excitation rate (at 
the final temperature - see below for a discussion of this 
point). The excitation time and the non-thermal velocity 
were used to calculate a free path travelled by the ion before 
it is excited. The final temperature of excitation was then 
computed using a mean temperature gradient. Using val- 
ues of the non-thermal velocity and the temperature gradi- 
ent appropriate to different regions of the atmosphere (quiet 
Sun, coronal holes, active regions), she found intensity en- 
hancements of up to a factor of 5. This was of the right 
order needed to explain the discrepancy bet ween o bserved 
and calculated intens ities r eported in Jordan (1975). 

Andretta et al. (2000) followed the same approach for 
the He II resonance line, but instead of a using single up- 
ward non-thermal velocity, they considered an ensemble of 
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plasma elements (clumps) in an isotropic turbulent velocity 
field with a Gaussian distribution having a standard devi- 
ation equal to the root-mean-squared non-thermal velocity 
yj{v\). Andretta et al. (2000) also used 7] = 8 x 10 4 K and 
assumed a constant pressure in the transition region. They 
calculated the enhancement factor over a parameter space 
of and dT /dh appropriate to conditions in the TR, and 
for values of the electron pressure corresponding to quiet 
Sun and active region conditions. Enhancement factors of 
up to about 5 were obtained for values of the parameters 
typical of current models and observations of the quiet Sun 



transition region (Pe/k - 



K, v /(4) > 10 km 



s , dT/dh > 10 K km ). Scaling arguments derived from 
their analysis show reasonable agreement with observations 
of the Hen 303. 8-A line compared with the Om 599. 6-A 
line. They concluded that the process could account for at 
least some of the enhancement required in the 303. 8-A line 
intensity. 

Our calculations of the effects of non-thermal motions 



follow the approaches of Jordan (1980) and Andretta et al. 
(2000), but with certain changes. The above studies used 
the 'coronal approximation,' assuming that in the resonance 
lines radiative decay is balanced by collisional excitation. 
Radiative transfer calculations suggest that th is is certainl y 
realistic in the case of the He II resonance line (Smith 2000), 
so this assumption is also made here. In treating the He I 
lines, advantage is taken of information from radiative trans- 
fer calculations made using full atmospheric models. 

In the coronal approximation the enhancement factor in 
the intensity of the 303. 8-A line is given approximately by 
the ratio of the collisional excitation rate at the final tem- 
perature to that at the initial temperature. The collisional 
excitation rate dj is given by the expression 



a 



t3 (T e ) =8.63 x 1Q~ 6 ^-T' 



/jVeexp^-— j s 

where <?; is the statistical weight of the lower level and flij is 
the collision strength, which is approximately constant over 
the temperature range of interest. The intensity enhance- 
ment ratio is therefore given approximately by 



(1) 
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Assuming the electron pressure to be constant in the tran- 
sition region, using equation (|l]), this becomes 
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The final temperature Tf depends on the velocity of the 
emitting plasma in the direction of the temperature gradi- 
ent. In order to generalize the calculation of the enhance- 
ment factor to an ensemble of plasma elements Andretta et 
al. (2000) convolved equation (^) with a normalized Gaus- 
sian distribution of turbulent velocities with standard devi- 
ation yj (v 2 ^): 



G(v) 



exp 



v 
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to give an intensity ratio of 



G{v)dv. 



(4) 
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2.2 New formulations 

Equation (|B|) was used in calculations of the enhancement 
factors in the present work, but the free path and final tem- 
perature Tf of each plasma element were computed by differ- 
ent methods to those used by Jordan (1980) and Andretta 
et al. (2000). These new techniques use the temperature gra- 
dient (now allowed to change over the path of the clump of 
plasma) derived from emission measure distributions found 
from observations of other transition region lines. Specific 
values of (wx)(Tj) suggested by observations were used. For 
the 584. 3-A and 537.0-A lines of He I the excitation rate 
ratio in equation was also modified. 

Both Jordan (1980) and Andretta et al. (2000) defined 
the lifetime of the ground state as r = (Cij(Tf))~ . This 
provides a lower limit on the lifetime and hence on the final 
temperature reached by the moving plasma. This may be ap- 
propriate if ions are carried rapidly to a higher temperature 
and remain there until excited. Here we use t = (Cij(Tj))~ , 
as we consider this timescale to better reflect the dependence 
of the enhancement process in question on the long excita- 
tion and ionization times of the helium ions at the equilib- 
rium temperatures of resonance line formation. It is these 
timescales, long with respect to those for other TR lines, 
which allow unresolved turbulent motions to carry helium 
ions to the higher temperatures where excitation occurs. 

7] is chosen to represent the 'normal' temperature of 
resonance line formation, in the absence of non-equilibrium 
processes. Tf is found by integrating the temperature gradi- 
ent over the path traversed in time r. The (varying) temper- 
ature gradient between Tj and Tj is derived from the mean 
intr insic emission measure distribution (EMD; see Section 
2.3). The intrinsic emission measure is defined as 



£m(0.3) = / N e N H dh (6) 

J Ah' 

where Ne and iVn are electron and hydrogen number densi- 
ties respectively and Ah' is defined as the interval of height 
in a plane-parallel atmosphere corresponding to A (log Te) = 
0.3 dex, the temperature interval over which most lines are 
typically formed. Defining electron and hydrogen pressures 
as Pe = NeTe and Pjj = N^Te, the emission measure may 
be re- written as an integral over the temperature interval in 
which a line is formed: 



Bm(0.3) 



s/2T' 



PePn dh 
dT e 



T 2 



dT e . 



(7) 



Te is the temperature at which the line contribution function 
peaks, and is the geometric mean of the limits of the integral, 
which represent a change in Te of a factor of 2. If Pe, Pu, 
and dh/dTe are each averaged over the temperature range 
of the integral, the emission measure becomes 



£m(0.3) = 



PePn dh 
V2T e dTe ' 



(8) 



The assumption of constant pressures over the region of line 
formation is a good one, and although the assumption of a 
a constant temperature gradient in the same region is less 
good, it is often used in constructing an initial model from 
the EMD (see e.g. Jordan & Brown 1981). 

If the mean EMD as a function of temperature is known 
from several lines, the path of material travelling up the tem- 
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perature gradient can be divided into smaller temperature 
intervals. Here we use intervals of 0.1 dex in log Te. In each 
of these intervals, the mean temperature gradient is calcu- 
lated from the EMD, at the mean value of log Te, using 
equation (|). In most of the temperature range of interest, 
hydrogen is fully ionized and Nn = 0.83Ae. This is a very 
good approximation above log Te = 4.9 (i.e. in the He II cal- 
culation), and is correct to better than 10 per cent down to 
log Te = 4.5 (the temperature range of the He I calculation). 
Thus, in plane-parallel geometry, taking Ph = 0.83Pe, the 
increment of height Ah corresponding to the temperature 
interval A (log Te) is given approximately by 



Ah = Em(0.3) 



V2T e 



AT e 



(9) 



where ATe spans the range log Te ± 0.05. The travel time 
required by a clump of material to traverse each Ah is de- 
termined by assuming a constant upward velocity v. The 
definition of Ah thus introduces a quantization of the (loga- 
rithmic) temperatures used in calculating the enhancement 
factor. For computational convenience this is carried over 
into the integration over velocities in equation (||). In this 
formulation, the integrand may only be evaluated correctly 
for values of Tc occurring at intervals of 0.1 dex in log Te 
from the initial temperature Tj (log Tc = log Tj + O.ln, for 
integer n). The values of v used as points in the integration 
correspond to velocities required to reach these Tc in a time 
equal to the lifetime of the ion's ground state at T:. The 
velocity mesh that results introduces some uncertainty into 
the calculation of the enhancement factors. 



2.3 Atmospheric parameters 

The EMDs used here to derive the variation of tempera- 
ture with height in the lower TR were determined from ob- 
served line fluxes and the appropriate atomic data using 
now-standard methods (see MJ99 and references therein). 
MJ99 derived separate EMDs for mean network and cell in- 
terior regions from observations of UV lines from the quiet 
Sun transition region made using the CDS (Coronal Diagnos- 
tic Spectrometer) and SUMER (Solar Ultraviolet Measure- 
ments of Emitted Radiation) instruments on board SOHO. 
Equation (^), used to determine the temperature gradi- 
ent from the emission measure, assumes a one-component 
plane-parall el atm osphere. In the network model proposed 
by Gabriel $97$ , open magnetic field lines are most ver- 
tical over the boundaries, but expand into the corona with 
a horizontal component in the cell interiors. Since we adopt 
a plane parallel atmosphere we use the network boundary 
EMD of MJ99 as the observed EMD. 

Of the lines observed by MJ99, those formed at tem- 
peratures below 10 5 K (except for He n) were observed with 
SUMER, and those formed at higher temperatures (and the 
helium lines) were observed with CDS. MJ99 found that the 
relative intensities of the He I lines observed with CDS and 
with SUMER in second order suggest that the SUMER lines 
might need to be redu ced in intensity by a factor of 1.5 if 
the Landi et al. ( 1997 ) CDS calibration were adopted. This 
would lead to the emission measure distribution below 10 
K, derived from the SUMER lines, being lowered relative to 
the higher temperature part derived from the CDS lines. 
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Figure 1. Network emission measure distributions. For log Te < 
5.3, the solid line shows EMD S, the dotted line shows EMD 
X. Theoretical EMDs are shown for log T e > 5.3, with (a) 
-Pe(log Te = 5.3) = 2 X 10 14 cm~ 3 K (dashed), and with (b) 
P e (log T e = 5.3) = 6 X 10 14 cm" 3 K (dot-dashed). 



This uncertainty inspired the use of a second EMD, de- 
rived from observations of £ Boo A, a G8 V star, made with 
the Goddard High Resolution Spectrogra ph (GHRS) instru- 
ment on the Hubble Space Telescope (Smith 200C). The 



EMD was scaled by a constant factor such that its minimum 
(at 1 - 2 x 10 5 K) fits the minimum of the network emission 
measure of MJ99. The resulting EMD is lower, by up to an 
order of magnitude, at temperatures below log Te = 4.8, 
but is very similar t o the solar EMD at higher tempera- 



tures. Jordan et al. (1987) found that below T e ~ 10 5 K 



the emission measure distributions of main sequence stars 
with hot coronae have approximately the same dependence 
on Te, which is supported by a direct linear correlation be- 
tween the C II and C IV emission line fluxes in such stars 



( |Rutten et al. 199l| ). Capelli et al. ( |1989[ ) found that the 
C II - C IV correlation also holds for solar observations at 
various spatial resolutions. 

The discrepancy between the SOHO observations and 
the earlier results may be due in part to the different spatial 
resolutions of the two instruments used (CDS and SUMER). A 
more serious problem is that SUMER could not observe the 
same location on the Sun in all of the lower TR lines, or 
the same location at the same time as in the CDS observa- 
tions (see Jordan et al. 2001). Thus there are fundamental 
uncertainties in the shape of the EMD derived from SOHO 
observations at temperatures below log Te — 4.9. The £ Boo 
A EMD was therefore used to explore the possible effects of 
this uncertainty on the present calculations. The solar EMD 
derived by MJ99 is hereafter referred to as EMD S, and that 
based on the £ Boo A EMD is referred to as EMD X; they 
are shown in Figure |l| 

The lines observed by MJ99 did not allow determination 
of the full EMD above log T e = 5.3. The EMDs used in the 
present calculations above this temperature are therefore 
based on energy balance arguments described below (see 
also Jordan 2000 and Philippides 1996). 

In the lower TR, below log Te — 5.0, some non-thermal 
source of heating is required to account for the rising emis- 
sion measure at low temperatures. In the upper TR, how- 
ever, the radiation losses are small compared with the con- 
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ductive flux from the corona, and may be balanced by energy 
deposited by the divergence of the conductive flux, or, if that 
flux is conserved, by a small amount of additional heating. 
The EMDs used here were determined assuming the former 
to be the case above log Te = 5.3. The equation expressing 
energy balance as a function of height in this region is then 

1 d(A(r)F c ) 



A(r) 



dr 



(10) 



where is the rate of radiative energy loss per unit volume, 
and F c is the classical conductive energy flux through an 
area A(r) normal to the radial co-ordinate r. The radiation 
loss term is given by 



e R = NeN K P l;id (Te) 



(11) 



where P rg A(Te) is the radiative power loss in all lines and 
continua, as a function of temperature in a plasma of known 
composition. P la A (Te) is calculated here using the power law 
f it to the radiative power loss curve adopted by Philippides 
([19960 : 



P rad (T e ) = 1.25 x 10" 
Fc is given by 



/ e erg cm s 



F C (T e 



-kT, 



5/2 



die 
dr 



(12) 



(13) 



where the Spitzer-Harm conduct ivity k ~ 1. 1 x 10 6 erg 



K 



-7/2 



at T e 



10 6 K flSpitzer 1956; ). That value 
is used here, as k is only 30% smaller at Te = 10 5 K. 



En- 



ergy transfer by mass motions of any kind (including the 
turbulent motion assumed in the enhancement factor cal- 
culations) is not considered in the energy balance equation 
used in the upper TR. 

Combining the energy balance equation ( |Io| ) with equa- 
tion (^) for the emission measure in terms of the temperature 
gradient gives an expression for the loga rithmic gradi ent of 
the emission measure with temperature (Jordan 2000): 

d log £m(0.3) 



d log Te 



3 dlogPe dlog A(T e ) 
2 dlogTe dlogTe 
2£m(0.3) 2 P rad (T e ) 



0.83/tPiTp 72 



(14) 



Assuming forms for P e (Te) and A(Te), this equation 
may be integrated numerically from given boundary con- 
ditions to give the emission measure as a function of tem- 
perature in that part of the atmosphere where the energy 
balance is described by equation (hoi) . In this case the elec- 
tron pressure was assumed to vary according to hydrostatic 
equilibrium (the use of a constant value of Pe in equation 
(^) in calculating the intensity enhancement is a good ap- 
proximation in the lower TR, owing to its small vertical ex- 
tent). A(Te) was taken only to vary with the expansion with 
height of a spherically symmetric atmosphere; assuming that 
emission comes from a uniform layer at each height, this is 
consistent with the use of the empirical EMDs derived in 
the plane-parallel approximation in the lower TR, because 
the small extent of that part of the atmosphere makes the 
two approaches approximately equivalent. 

The assumptions about A(Te) represent a significant 
simplification. The EMD derived below log Te = 5.3 rep- 
resents an average over any unresolved structure within 
the network boundaries, while the EMD computed using 



equation (^) is the 'intrinsic' value; they are unlikely to 
match at log Te = 5.3. If the calculated emission measure at 
log Te = 5.3 is greater than the value derived from spatially 
averaged observations, this implies that only some fraction 
of the observed region is occupied by emitting material. 

With the assumed forms for P e (Te) and A(Te), the 
EMD for the upper TR was determined using equation (|l4j) , 
choosing values of the emission measure, electron pressure, 
and coronal temperature Tq at the top of the region. Tq 
was chosen such that log Tq = 6.20, as CDS observations 
suggest that there i s little material at hi gher temperatures 
in the quiet corona ( Landini & Landi 199£ ) . For a given coro- 



nal pressure and temperature, the coronal emission measure 
cannot exceed a certain value if the RHS of equation ( |l4| ) 
is to remain positive above log Te = 5.3, as is required by 
the observed gradient of the EMD. There is hence a 'criti- 
cal solution,' with the largest possible value of the coronal 
emission measure which reproduces the observed minimum 
in the EMD at log T e ~ 5.3. If Em(T c ) is larger than this 
value, conduction from the corona deposits more energy in 
the upper TR than can be lost by radiation. 

The mean pressure in the boundary regions used in de- 
riving EMD S below log T e = 5.3 was found by MJ99 to be 
log Pe = 14.27, with a standard deviation of ±0.35. Using a 
coronal pressure that reproduces this pressure (~ 2 x 10 14 
cm -3 K) at log Te = 5.3 in the critical solution gives EMD 
(a) shown in Figure |l| This is lower than the observed emis- 
sion measure at log Te = 5.3, suggesting that a slightly 
higher pressure is needed for a continuous EMD through 
this temperature. EMD (a) also predicts substantially less 
Mg IX and Mg x emission than observed by MJ99 (when 
the photospheric abundance of Mg is used in both P I&c \(Te) 
and in calculating the line intensities). There are problems 
with such a comparison, however, as the intensities used by 
MJ99 correspond to areas defined as network boundaries 
in the CDS transition region lines, and the details of the 
connection between the TR network and the corona are un- 
known. Open field regions would be expected to expand in 
the corona, but the network boundaries may also contain 
closed field emitting regions unconnected to the corona. 

Figure | also shows EMD (b) , a second upper TR EMD 
calculated to give a higher transition region pressure of Pe = 
6 x 10 14 cm" 3 K, which corresponds to the largest pressures 
found by MJ99 (in cell interiors) and is similar to values 
used in earlier models of the atmosphere (e.g. Vernazza et 
al. 1981). It provides a much better match to the Mg IX and 
Mg x data, although, as stated above, these do not provide a 
stringent constraint. EMD (b) predicts an intrinsic emission 
measure at log Te = 5.3 that is larger than the apparent 
(plane-parallel) emission measure. This discontinuity can be 
explained in terms of an area filling factor below log Te = 5.3 
(see Sim & Jordan 2001). For a filling factor A^/A Q u s , the 
intrinsic emission measure below log Te = 5.3 is 



Em: = 



obs 
A: 



Em&pp 



(15) 



where -Emapp is shown in Figure In the theoretically de- 
rived upper TR the filling factor is assumed to be 1.0. Scal- 
ing the lower TR EMD to match the upper TR EMD (b) at 
log Te = 5.3 suggests a filling factor in the lower TR of 0.16 
(within the area of ~ 3.4 arcsec 2 resolved by CDS). The con- 
tinuous EMD then describes the restricted areas where the 
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filling factor is 1.0. The lines used to find the lower TR EMD 
are not sensitive to Pe, so a similar matching procedure may 
be considered for any lower pressure above about 2.7 x 10 14 
cm -3 K, for which the theoretical upper EMD would match 
the observed emission measure at log Te = 5.3, with a filling 
factor of 1.0. Equations Q and (jl5]) show that the temper- 
ature gradient derived from Em(Te) scales as P^A^/A ^ S , 
but the value of A^/A^g (i- e - the factor by which the theo- 
retical and observed EMDs differ at log Te = 5.3) required 
for a continuous EMD through log Te = 5.3 scales as PjT 2 . 
If the filling factor in the lower TR is regarded as a free 
parameter, the distribution of dTe/dh is the same for all 
continuous EMDs in which a lower TR EMD is scaled to 
match an upper TR EMD, and it is independent of the as- 
sumed pressure. 

The distributions of dTe/dh derived in this manner for 
the cases of lower TR EMDs S and X were used in the cal- 
culations of the helium line enhancement factors. In order 
to investigate the effect of the possible solutions suggested 
by EMDs (a) and (b), the calculations were performed for 
electron pressures representing these limits, Pe = 2 x 10 14 
cm" 3 K and Pe = 6 x 10 14 cm" 3 K (the lower pressure is 
not technically consistent with the filling factor argument, 
but, considering other approximations involved in the cal- 
culations, provides a useful lower limit). Because the tem- 
perature gradients used in the enhancement factor calcula- 
tions are independent of the assumed pressure, the pressure 
adopted influences only the excitation time. The excitation 
time scales as P^ 1 ; hence the velocity required to reach a 
given Tf, v(Tc), scales as Pe- A higher pressure will there- 
fore lead to larger v(Tt), causing contributions from higher 
temperature regions to be suppressed by the Gaussian term, 
G(v), in the enhancement factor integral. 

The values of (v^) used in the integral in equation ^ 
were determined from empirical relations fitted to observa- 
tional data. From early rocket experiments (e.g. Berger et al. 
1970; Boland et al. 1973) to recent results from SOHO (e.g. 
Chae et al. 1998; Peter 2001), solar observations of optically 
thin TR lines have shown widths in excess of those due to 
thermal motion at the temperatures at which the lines are 
thought to form. The non-thermal component of the broad- 
ening could be due to unresolved upflows and downflows, to 
acoustic or MHD waves, or to small scale turbulence. We 
interpret the non-thermal widths in the TR as evidence of 
turbulent motions of unspecified origin. Since the line pro- 
files are approximately Gaussian, assuming a Gaussian dis- 
tribution of non-thermal velocities allows the most probable 
turbulent velocity £ to be extracted from the line width us- 
ing the expression 



AA 
A 



2fcT; 



ion _|_ ^2 



1/2 



(16) 



where AA is the full width at half maximum, A is the wave- 
length of the line, Pj on is the ion temperature of line for- 
mation, and mj on the mass of the emitting ion. The mean 
squared non-thermal velocity {v^} is given by the relation 
(4> = (3/2)£ 2 . 

Up to log Te ~ 5 .3, both recent SOHO measurements of 



£ (Chae et al. 1998) and earlier measurements (see Jordan 



2 ~ rrl/2 which 



conserved through the transition region to be dissipated at 
greater heights. The values of (vy)(T;) used in the present 
study are found using 



/ T \ 1/4 

^H^k) kms 



(17) 



The value of £ 2 and hence (v^) used to characterize the 
velocity field in each calculation was taken to be the value 
given by e quatio n ( p^ ) at the initial temperature Tj. 

Peter (2001) has found that some transition region lines 
(forming in the temperature range 4.6 < log Te < 5.8) can 
be fitted more precisely with two Gaussian components. The 
widths of the core components are similar to those found by 
Chae et al. (1998); the widths of the broad components are 
signific antly larger, but are also consistent with £ 2 oc T ( 



1/2 



Peter (2001) interprets the broad components as evidence 
of the passage of magnetoacoustic waves in coronal fun- 
nels in the network. It is not known whether the (optically 
thick) helium lines have this broad component. If larger tur- 
bulent velocities were assumed in the present calculations, 
contributions from higher temperatures would become more 
important and the enhancement factors derived would be 
larger. The temperature dependence of £ 2 would mean that 
the enhancement factors for He II would be increased by a 
relatively greater amount than those for He I. 



3 ENHANCEMENT FACTOR CALCULATIONS 
3.1 The He II 303. 8-A line 

Radiat ive transfer calculations with a 39 level model helium 
atom ( [Smith 2000 ) show that the contribution function to 
the intensity of the 303. 8-A line peaks at log Te — 4.9, and 
this value was used for log in the integration of equation 
(^) for the He II line. The same peak formation tempera- 
ture was found by MJ 99 usi ng the ionization equilibrium of 

assuming purely collisional ex- 



Arnaud & Rothenflug (198- 



citation, which is consistent with the model calculations. At 
log T- = 4.9, y/{v^)(T{) = 24.3 km s _1 . The collisional ex- 
citation rate Cij(Te) for the He II resonance transition was 
calculated using equation ([[]), taking Pe as constant and us- 
i ng th e collision strength fiy(Te) given by Aggarwal et al. 
(|199§). This rate was also used to calculate the lifetime of 
the ground state, r = Cy 1 . The He II excitation time at this 
temperature is 5.69 s, compared with a time of 0.06 s for the 
O m 5 99. 6- A line, using the collision strengths of Aggarwal 




and 3 x 10" 3 s for the C IV 1548- A lines (Mariska 



1991) are fitted quite well by the relation £ oc T e 
would be expected if a wave flux from lower regions were 



formed at slightly higher temperatures of log Te — 5.0 
(these times assume Pe = 6 x 10 14 cm -3 K). The differences 
imply that other lines formed at similar temperatures to the 
He II 303. 8-A line would be relatively unaffected by turbu- 
lent transport. 

In order to facilitate a discussion of some of the inter- 
esting features and limitations of the enhancement factor 
calculations, details of the calculations are presented in Ta- 
ble [l] and Figures u and |. Table shows the heights above T: 
corresponding to the final temperatures Tf determined from 
the EMDs using equation (|^). Figure [| shows the excitation 
rate ratio as a function of log Te, and Figure ^ illustrates the 
form of the integrand in equation Q , showing the contribu- 
tion to the total intensity enhancement of plasma elements 
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Figure 2. The excitation rate ratio Cij(Tf)/Cij(T:) as a function 
of (logarithmic) temperature for the He II 303. 8-A line, as used 
in calculations of the enhancement factors. Log Tj = 4.9. 



Table 1. Final temperatures and heights above log = 4.9 
reached by plasma elements with velocity v(Tf), computed us- 
ing EMDs S and X. An assumed pressure of Pe = 2 X 10 14 cm -3 
K is denoted by (a), P e = 6 X 10 14 cm -3 K is denoted by (b). 
Values of Tf < T; correspond to downward-moving elements. 



log T f 



h(T t ) - h(3i) 
(km) 



v(T { ) 
(km s" 1 ) 





S 


X 


S(a) 


X(a) 


S(b) 


X(b) 


4.6 


-84.4 


-44.8 


-4.94 


-2.62 


-14.8 


-7.87 


4.7 


-38.2 


-33.7 


-2.24 


-1.97 


-6.71 


-5.92 


4.8 


-13.0 


-17.3 


-0.762 


-1.01 


-2.28 


-3.04 


5.0 


11.6 


12.3 


0.680 


0.721 


2.04 


2.16 


5.1 


28.0 


25.3 


1.64 


1.48 


4.92 


4.45 


5.2 


51.2 


45.9 


3.00 


2.69 


9.00 


8.07 


5.3 


83.9 


78.6 


4.92 


4.60 


14.7 


13.8 


5.4 


110 


135 


8.20 


7.91 


24.6 


23.7 


5.5 


262 


257 


15.3 


15.1 


46.0 


45.2 


5.6 


518 


513 


30.3 


30.1 


91.0 


90.2 


5.7 


1070 


1060 


62.4 


62.1 


187 


186 


5.8 


2190 


2180 


128 


128 


384 


383 



with different velocities. The velocities required to reach Tf 
are listed in Table [lj The quantization of log Te and h pro- 
duces a fine mesh of velocities at and below the peak in 
the excitation rate ratio, but a coarser grid at higher tem- 
peratures. This results in some over-estimation of the total 
He II enhancement factor, particularly in the contribution 
from the material with v > 30 km s . Rough calculations 
suggest this over-estimation could be by about 10-20 per 
cent. 

The calculations made assuming Pe = 6 x 10 14 cm~ 3 
K show the peak contribution to the intensity occurring for 
upward velocities of about 15 km s _1 , but with a signifi- 
cant contribution from faster-moving material. The contri- 
bution of material moving downwards (v < 0) is seen to be 
negligible. The peak contribution to the intensity occurs at 
log T e = 5.3, at a height above Tj of about 80 km. EMDs 
S and X lead to similar temperature profiles, resulting in 
similar enhancement factors being derived for each, 5.3 for 
EMD S and 5.5 for EMD X. 
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Figure 3. The integrand in equation (|5j) for the intensity en- 
hancement, evaluated for the He II resonance line, using EMD S 



(solid) and EMD X (dotted) with P e 
and P e = 6 X 10 14 cm" 3 K (bottom). 



2 X 10 1 



" 3 K (top) 



Assuming a lower pressure with these EMDs prod uces 
larger enhancement factors, as predicted in Section 2.3. The 
total enhancement factors computed for Pe = 2 x 10 cm -3 
K are 6.9 for EMD S and 7.0 for EMD X. Figure | and Table 
[jjshow that the peak contribution to the intensity occurs at a 
greater height and temperature than at the higher pressure, 
as the velocities required to reach given Tf are decreased 
(because the excitation time is longer). 



3.2 The He I 584. 3-A and 537.0-A lines 

Similar calculations were performed for the first two lines 
of the He I resonance series, although here the situation is 
more complicated, as the coronal approximation no longer 
applies. First, the He I resonance lines form at greater op- 
tical depths than the He II line, and full radiative transfer 
modelling is required to compute their intensities properly. 
Even neglecting transfer of photons within the line, direct 
collisional excitation is in general not the dom inant excita- 
tion process in radiative transfer calculations ( Hearn 196£ ; 



Andretta & Jones 1997; Smith 2000). Those calculations re- 



veal that important contributions come from radiative tran- 
sitions from levels other than the ground state, which are 
themselves populated by a combination of collisional excita- 
tion and radiative recombination. Neither the excitation rate 
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nor the mean lifetime of the He I ground state can be accu- 
rately evaluated using the equivalent expression to equation 
(^) used in the He II calculation. In order to investigate 
the possible effect of non-thermal motions on the 584. 3-A 
(ls2p 1 P - Is 2 1 S) and 537.0-A (ls3p 1 P - Is 2 1 S) lines in 
a first approximation, corrections were made to the expres- 
sions used to compute these quantities. Radiative transfer 
effects were not considered. 



Radiative transfer calculations ( smith 200C ) reveal that 
significant contributions to the 2 1 P level population come 
from allowed radiative transitions from other singlet levels, 
especially 2 1 S, 3 1 5', and 3 1 D. These levels are collisionally 
excited from the ground at high temperatures but largely 
populated by radiative recombination at low temperatures. 
The collisional contribution to the population of the 2 1 P 
level through these channels was approximated in the calcu- 
lation of the velocity redistribution enhancement factor by 
writing the total excitation rate CL(Te) as 



C ii (T e ) = C ij {Te) + J2 



A 



kj 



■C ik (T e 



(18) 



where CVj(Te) is the direct collisional excitation rate, and 
the second term is the sum of the collisional excitation 
rates of the other contributing levels, multiplied by radia- 
tive branching ratios. CV;(Te) and Cifc(Te) were evaluated 
using expression ([!]), with collision strengths for Te = 3x 10 4 
K (the collision strengths were assumed to be constant at 
higher temperatures, as their variation at higher Te w as of- 
ten unavailable) provi ded b y Lanzafame et al. (1993) and 
Sawey & Berringt on ( 1993 ) . The Einstein A- values were 
taken from Drake (1996). 

The excitation rate ratio was calculated using log Xj = 
4.5. This temperature was found by MJ99 to represent the 
peak of fo rmation of th e 584. 3-A line. Radiative transfer cal- 
culations (Smith 2000) suggest a slightly lower temperature 
log Tj ~ 4.45 due to a different ionization balance, but as 
comparisons are being made here with MJ99, the former fig- 
ure was kept. Using a lower T; would have relatively little 
effect on the final temperatures responsible for the velocity 
redistribution intensity enhancement, as at log Tj ~ 4.45, ex- 
citation times for the He I resonance line are so long (about 
40 s; Smith 2000) that ions in the ground state are very 
likely to survive at least to log Tj = 4.5, and their subse- 
quent behaviour to be determined by higher temperatures 
(as for He II, peak enhancement occurs where log Te > 5.0). 
The absolute value of the intensity enhancement is, how- 
ever, more significantly affected by the choice of Ti, through 
the excitation rate ratio. The total collisional excitation rate 
calculated using equation ( |l8| ) is a factor of two smaller at 
log Tj = 4.45 than at log T; = 4.5, so choosing a lower Tj 
would be expected to produce a larger enhancement. 

The expression used for the excitation rate is a less ex- 
act approximation to the true rate than in the case of He II, 
neglecting as it does the effects of recombination, which is 
proba bly important in the formation of th e He I resonance 
lines (Andretta & Jones 1997; Smith 2000). Recombination 



is most significant below the peak temperature of line for- 
mation. At log Te > 4.5, the temperature range considered 
here, the expression used represents about 75 per cent of 
the total excitation rate computed in radiative transfer cal- 
culations. The fraction does not change significantly with 




Figure 4. The excitation rate ratio ■ (Tf)/C? ■ (Tj) as a func- 
tion of (logarithmic) temperature for the He I 584. 3-A (solid) and 
537.0-A (dotted) lines, as used in calculations of the enhancement 
factors. Log Tj = 4.5. 



temperature in this range. As the enhancement factor due 
to non-thermal transport of He depends on relative rates at 
Tj- and Tj, this approximation should have little effect on the 
effect on the calculated values of the enhancement factors. 

Whereas in He II the collisional excitation of the res- 
onance line completely dominates the depopulation of the 
ground state, in He I there are competing processes which 
act to shorten the lifetime of the ground state, some of which 
do not result in excitation to the 2 1 P level. Collisional exci- 
tation rates to other bound levels are significant, as consid- 
ered above, and collisional ionization becomes important at 
high temperatures. Hence the probability that a given He I 
ion survives in the ground state until reaching Tf is smaller 
than would be determined by simply taking the reciprocal of 
C'ij (equation (^)). In order to take this shortening of the 
effective excitation time into account in a first approxima- 
tion, for He I t(Tj) was calculated by removing the radiative 
branching ratios from the expression for C'g and adding the 
collisional ionization rate (Mihalas & Stone 1965) before tak- 



ing the reciprocal to give a mean lifetime (no longer strictly 
an 'excitation time'). 

Using these new prescriptions for determining C^JT:) 
and r(Tj), the intensity enhancement integrals for the 584.3- 
A line were performed as for the He II resonance line, with 

\/i v T)( T i) = 19A km s_1 and r ( T i) = 10 - 5 s ( this time 
may be compared with the excitation times, at the same 

temperature, of the C II lines near 1335 A of about 0.05 
s - Smith 2000). Details of the integrals are presented in 
Figures ^| and ^ and Table ^. These show that assuming Pe = 
6 x 10 14 cm -3 K with EMD S predicts the peak contribution 
to the intensity to come from material with v ~ 20 km s _1 , 
reaching log Te ~ 5.1 at a height of about 200 km above 
Tj. Significant contributions are again seen from material 
reaching higher temperatures. Taking Pe = 2 x 10 14 cm -3 
K with EMD S, the peak contribution comes from slower- 
moving material (v < 10 km s _1 ), but because the excitation 
time is longer at the lower pressure, this material reaches 
greater heights and is excited at higher temperatures. The 
total enhancement factor is correspondingly increased: 15.3 
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Table 2. Final temperatures and heights above log T- = 4.5 
reached by plasma elements with velocity v(T^), computed us- 
ing EMDs S and X. An assumed pressure of Pe = 2 X 10 14 cm -3 
K is denoted by (a), P e = 6 X 10 14 cm -3 K is denoted by (b). 
Values of Tf < T; correspond to downward-moving elements. 



log T f h{T { ) - h(Ti) v(T f ) 

(km) (km s _1 ) 





S 


X 


S(a) 


X(a) 


S(b) 


X(b) 


4.2 


-181.5 


-30.3 


-5.77 


-0.963 


-17.3 


-2.89 


4.3 


-99.3 


-13.1 


-3.16 


-0.416 


-9.47 


-1.25 


4.1 


-58.1 


-5.95 


-1.85 


-0.189 


-5.54 


-0.567 


4.6 


92.4 


7.68 


2.94 


0.244 


8.81 


0.732 


4.7 


139 


18.8 


4.42 


0.597 


13.3 


1.79 


4.8 


164 


35.2 


5.21 


1.12 


15.6 


3.36 


4.9 


177 


52.5 


5.62 


1.67 


16.9 


5.00 


5.0 


188 


64.8 


5.97 


2.06 


17.9 


6.18 


5.1 


205 


77.8 


6.51 


2.47 


19.5 


7.42 


5.2 


228 


98.4 


7.24 


3.13 


21.7 


9.38 


5.3 


261 


131 


8.29 


4.16 


24.9 


12.5 


5.4 


317 


187 


10.1 


5.94 


30.2 


17.8 


5.5 


439 


309 


13.9 


9.82 


41.8 


29.5 


5.6 


695 


565 


22.1 


18.0 


66.3 


53.9 


5.7 


1240 


1110 


39.5 


35.3 


118 


106 


5.8 


2360 


2230 


75.1 


70.9 


225 


213 



for P e = 2 x 10 14 cm" 3 K compared with 11.2 for P e 
f i x 10 14 cm" 3 K. This can be seen clearly in Figure 

For both values assumed for the pressure, the total en- 
hancement factors computed using EMD X are larger than 
for EMD S. A factor of 17.6 is found for P e = 2 x 10 14 cm" 3 
K and a factor of 18.7 is found for Pe = 6 x 10 14 cm" 3 K. 
In both cases the peak in the intensity contribution occurs 
at log Te = 5.2 at a height of about 100 km above Tj. The 
velocities required to reach this If are smaller than in the 
calculations using EMD S, as the temperature gradients de- 
rived from EMD X in the lower TR are larger. The Gaussian 
weighting leads to the larger enhancements found with EMD 
X. The factors computed for He II using the two EMDs are 
more similar because the (lower TR) EMDs themselves are 
similar above the initial temperature (log 7] = 4.9). 

Although the use of log Tj = 4.5 for the He I resonance 
line may lead to under-estimation of the enhancement fac- 
tors if the line forms in equilibrium at a slightly lower tem- 
perature, other approximations used in the calculations tend 
to lead to over-estimation. In comparing with the work of 
MJ99, who assumed line formation at log Te = 4.5, the fig- 
ures derived here should therefore be regarded as upper lim- 
its on the possible enhancement factors. This can be seen in 
the fact that the enhancement of the 584. 3- A line is found to 
increase with pressure in calculations using EMD X, show- 
ing that the uncertainties introduced by the large velocity 
steps in the tail of the integrand are significant in this case. 

MJ99 did not compare the observed intensity of the He I 
537. 0-A line directly with the EMD derived from 'normal' 
transition region lines, but radiative transfer calculations 
under-produce the observed intensity of this line by factors 
of the same order as in the 584. 3- A line (Smith 2000). En- 
hancement factors were therefore computed for the 537. 0-A 
line, focusing on a comparison with the 584. 3-A line. Indi- 
rect collisional excitation (mainly through the 2 1 S, 4 1 S 
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Figure 5. The integrand in equation (|5j) for the intensity en- 
hancement, evaluated for the He I resonance line, using EMD S 
(solid) and EMD X (dotted) with P e = 2 X 10 14 cm" 3 K (top) 
and P e = 6 x 10 14 cm" 3 K (bottom). 

and 4 1 D levels) was again included in the total excita- 
tion rates calculated for T- (again log T- = 4.5) and Tf and 
in the calculation of the mean ground state lifetime. The 
537. 0-A line excitation ratio shows similar variation with 
velocity (and Tf) to that of the 584. 3-A line, but owing to 
the higher excitation energy of the 3 1 P level and the in- 
termediate levels important in the indirect excitation of the 
537. 0-A line, the excitation rate ratio peaks at a slightly 
higher temperature and falls off more slowly (relative to its 
absolute magnitude) than for the 584. 3-A line. Very simi- 
lar calculations to those performed for the 584. 3-A line were 
carried out. The integrated intensity enhancement factors 
computed for the 537.0-A line are 14 (EMD S) and 27 (EMD 
X) for P e = 6 x 10 14 cm" 3 K, and 22 (S) and 26 (X) for 
Pe = 2 x 10 14 cm -3 K. The trends exhibited are the same 
as those seen in the calculations for the 584. 3-A line. The 
ratio of the enhancements computed for the two He I lines 
is discussed in Section 0. 



3.3 Non-Maxwellian effects 

Many studies of the electron velocity distribution function 
(EVDF) in the solar transition region have suggested that 
significant departures from a Maxwellian distribution may 
occur (e.g. Shoub 1983; Ljepojevic & Burgess 1990; Vinas, 
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Wong & Klimas 2000). This could be due to streaming of 
fast electrons from the high TR and corona down the steep 
temperature gradient or due to accelerating processes in the 
chromosphere or transition region. A form of the EVDF ap- 
propriate to the former process has been used an extensive 
investigation of the effects of such a distribution on the he- 
lium spectrum (Paper II). Here we briefly consider the ef- 
fects of that non-Maxwellian EVDF on the enhancements 
predicted by the turbulent transport calculations. 

An EVDF that is Maxwellian below a c ut off velocity 
of 3.5 times the thermal velocity, yj2k Te/m, with a power 
law decline with v~ ~ 31/7 flShoub 1982| ) at higher velocities 
was considered. This represents a particularly large depar- 
ture from the Maxwellian form, and was chosen in order 



Table 3. Upper limits on intensity enhancement factors for the 
helium lines calculated for EMDs S and X. An assumed pressure 
of P e = 2 X 10 14 cm" 3 K is denoted by (a), P e = 6 X 10 14 cm" 3 
K is denoted by (b). 



to p rod uce the maximum effect on collision rates (Smith 
200c|)~K 

quantitative investigation of the effects of such a 
distribution on collision rates (and hence on ground state 
mean lifetimes) was performed by comparing excitation and 
ionization rates computed in radiative transfer calculations 
with and without the non-Maxwellian tail. 

This comparison showed that the power law tail in the 
non-Maxwellian distribution function does not have a sig- 
nificant effect on the collision rates at temperatures much 
above the normal peak temperature of formation (Tj). At 
higher temperatures the excitation energies of the helium 
resonance lines correspond to the energies of electrons in the 
Maxwellian part of the local distribution, and these electrons 
dominate the excitation rates. At and below Tj, collisional 
excitation and ionization occurs predominantly by electrons 
in the suprathermal tail of the distribution, and so rates are 
significantly enhanced compared with the Maxwellian case. 
In the enhancement process investigated here, as excitation 
occurs mainly at temperatures much higher than Tj, the 
power law tail has little effect except at Tj. The principal 
effect is to reduce the mean lifetime of the He I and He II 
ground states at their respective Tj. For He n this effect is 
minimal, reducing r(Tj) by 10 per cent, but for He I the 
reduction is more significant, being by a factor of about 7. 

The increase in CV,(Tj) in equation (^) would also tend 
to reduce the computed enhancement factor, simply because 
emission at T\ would be increased; the effect would not be 
related to the efficiency of turbulent transport. The change 
in Cy(Tj) is therefore ignored in assessing the effects of the 
non-Maxwellian EVDF on the enhancement factor. 

The decrease of r(Tj) has negligible effect on the en- 
hancements computed for the He II line, but noticeable ef- 
fects on those computed for the He I lines. Excitation of the 
584.3-A and 537.0-A lines occurs mainly at the same tem- 
peratures and heights as in the calculations with Maxwellian 
excitation, but the velocities required by plasma elements to 
reach those heights are greater. There is, therefore, some re- 
duction of the total enhancement (the same effect occurs for 
He II, but to a very much smaller extent). The reduction is 
relatively greater in the case of EMD S (the enhancement 
factor is smaller by about 30 per cent) than for EMD X, as 
the temperature gradient above log Te = 4.5 is smaller in 
the former case. The velocity required to reach a particular 
temperature is increased for EMD S by a greater amount 
relative to U (v^,) (Tj) than for EMD X, reducing the Gaus- 
sian weighting G(v) to a greater extent in calculating the 
EMD S enhancement factor. 

In the 537.0-A line, the effect of the smaller G(v) at high 



Ion 


Wavelength 


Enhancement factor 




(A) 


S(a) 


X(a) 


S(b) X(b) 


He II 


303.8 


6.9 


7.0 


5.3 5.5 


He I 


584.3 


15.3 


17.6 


11.2 18.7 


He I 


537.0 


22.4 


26.0 


14.1 26.5 



temperatures leads to a reduction in the temperature of the 
peak intensity enhancement so that it occurs in the same 
region as for the 584.3-A line. This would tend to reduce any 
variation in the ratio of intensities of the two lines compared 
with the case in which Maxwellian excitation is assumed. 
The enhancement of the 537.0-A line is also decreased, by 
almost a factor of two in calculations with EMD S, but by 
a much smaller amount in the case of EMD X. 



4 DISCUSSION AND CONCLUSIONS 

The intensity enhancement factors calculated for the helium 
resonance lines in the presence of turbulent transport are 
given in Table |^. As discussed above, these should be re- 
garded as upper limits. In the case of He II the over-estimates 
should be relatively small, and an intensity enhancement of 
approximately a factor of 5 is plausible, given the parame- 
ters assumed to describe the solar atmosphere. This is sim- 
ilar to the factors predicted by Jordan (1980) and Andretta 
et al. (2000) using a less satisfactory expression for the ex- 
citation time. An enhancement factor of 5 cannot account 
entirely for the enhancement of at least a factor of 13 ap- 
parently required in the network according to the analysis 
of MJ99. However, using the current second order correc- 
tion of sim25 for the NIS 2 waveban d of C DS instead of the 
factor of 55 derived by Landi et al. (1997), the required en- 



hancement is only a factor of 6 (as found by Jordan 1975), 
which is within the range calculated here. Radiative trans- 



fer ca lculations (Smith 2000) using the VAL ( Vernazza et al 



1981) models of the quiet solar atmosphere under-produce 
the 303. 8- A line intensity by factors of 3-4, which the present 
work shows could be explained entirely by the enhancement 
mechanism suggested here. 

The enhancement factors derived for the He I 584.3-A 
line are also upper limits, but show that enhancement factors 
of order 10 are plausible. Thus non-thermal transport of 
He I could produce the enhancement of the resonance line 
intensity required to account for the results of MJ99. 

Calculations of the 584.3-A line intensity using the VAL 
C (average quiet Sun) and VAL D (average net work) models 
produce fa i rly good ma tches to observations (Andretta & 
Jones 1997; Smith 2000j ) , but produce too much emission in 
other low transition region lines, owing to the presence of the 
temperature plateau at Te ~ 2.5 x 10 4 K. However, models 
of the atmosphere in which the plateau is removed do not 
produce high enough intensities in the He I line to match the 
observations, but by smaller factors than required by MJ99 
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( Andretta fc Jones 1997 1_ Smith 200C ). These discrepancies or even O m (5.05), O IV (5.25) or O v (5.4) ( Bruecknei 



could a ^s^Jac^racpkaincd^b^v^hc^fc^^ 



Bartuo 19741; |Withbroe 1976) ; |Gallagher et al. 1998 



Pat- 



Sir jiilar arguments hold for the He I 537.0-A line, for sourakos et al. 1999[ ). Patsourakos et al. (1999) found that 



which an enhancement of a factor of up to about 25 is pre- 
dicted by the calculations presented here. It is interesting 
that greater enhancement factors are predicted for the 537.0- 
A line than for the 584. 3-A line, as radiative transfer calcula- 



tion s us ing the VAL models (Andretta & Jones 1997; Smith 
200c|)~s}iow the predicted ratio 1(537.0 A)/J(584.3 A) to be 
smaller than observed. MJ99 found a ratio of 0.116 ± 0.015 
in the network (using the Landi et al. 1997 calibration of 
CDS; 0.105 ±0.014 using the Brekke et al. 2000 calibration), 
but calculations using the VAL models give values of about 
0.08. Enhancements of the two lines in the ratio derived 
here would help resolve this disagreement. Radiative trans- 
fer calculations using a model atmosphere based on EMD 
S produce a value of 0.117 (Paper II); the enhancements 
predicted here would therefore worsen agreement with ob- 
servations. Nevertheless, the enhancement factors calculated 
here give an intensity ratio much closer to that observed than 
do radiative transfer calculations of the effects of non-local 
suprathermal electrons which produce enhancement factors 
of the same order (Paper II). 

In the formulation used here, the efficiency of the non- 
thermal transport process increases with decreasing pres- 
sure. For example, the predicted 303. 8-A and 584. 3-A line 
intensity enhancements are increased by about 30 per cent 
when the pressure is reduced by a factor of 3. This would 
appear to conflict with observations of coronal holes, where 
electron pressures are lower than in the quiet Sun by fa ctors 
of 2 - 3 ( |Munro fc Withbroe 1972; ; Jordan et al. 200l[ ), but 
where helium line intensities (and required enhancements) 
are smaller. However, the energy balance adopted here in the 
upper TR (equation ([h])) takes no account of the fast solar 
wind and the term in A(r) does not allow for super-radial 
expansion. For non-thermal transport calculations appropri- 
ate to coronal holes, energy balance calculations specifically 
for coronal holes are required; these could be combined with 
observed EMDs in a similar manner to that described here. 
Further observations of the helium lines in coronal holes on 
the disk would be useful (to avoid the problems of limb dark- 
ening/brightening seen in polar coronal holes). It would also 
be interesting to investigate non-thermal transport in active 
regions (for which Andretta et al. 2000 predict smaller en- 
hancements), where electron pressures are higher and tur- 
bulent velocities lower than in the quiet Sun. 

Turbulent transport could help explain the spatial vari- 
ations of the intensities of the helium lines with respect to 
each other and to other TR lines in the quiet Sun. The in- 
creases in the heights of peak line formation predicted here 
are not large enough to be resolved directly at the limb, but 
the effects of the increased temperatures of emission may 
perhaps be seen in observations of the network, the appear- 
ance of which changes between the low and high TR. The 
turbulent motions of the largely ionized gas would be influ- 
enced by the direction of the network field, so that material 
moving to greater heights would trace the expansion of the 
network (see e.g. Gabriel 1976). The width and the contrast 
of the network observed in the He I and Hell lines seem 
to resemble that seen in high TR lines like those of Ne VI 
(log T e = 5.6) and Ne VII (log T e = 5.7), rather than that 
seen in lower TR lines of C II (log T e = 4.4), C ill (4.9), 



Gabriel's (1976) model predicts a network width which in- 
creases with temperature, reaching the width observed in 
helium at log Te — 5.75. The present calculations predict 
He i emission to peak at log Te — 5.2, and He n emission 
to peak at log Te ~ 5.3, but with significant contributions 
in each case from temperatures up to log Te ~ 5.6. There is 
little contribution to intensity above log Te ~ 5.7, but the 
calculations assume different pressures and temperature gra- 
dients (and overall geometry) to those in the Gabriel (1976) 
model, making quantitative comparisons difficult. Qualita- 
tively, however, turbulent transport in an expanding net- 
work could explain why the pattern of helium emission re- 
sembles that seen in lines formed in equilibrium at higher 
temperatures and heights. In a similar geometry, suprather- 
mal electrons streaming down the field lines from the upper 
TR would be expected to be funnelled into the network in 
the low TR, which would concentrate any enhanced helium 
emission in the centre of the network, which would be less 
consistent with the observed width of the network. 

Turbulent transport of helium in an expanding net- 
work is also qualitatively consistent with observed vari- 
ations in the ratio of the intensities of the He I lines, 
7(537.0 A)//(584.3 A). The excitation rate ratio for the 
537.0-A line peaks at a slightly higher temperature and falls 
off more slowly at higher temperatures than for the 584. 3-A 
line. This suggests that the network could appear wider in 
the 537.0-A line than in the 584. 3-A line, which is broadly 
consistent with rastered images in the two lines (MJ99), al- 
though further observations of the two lines would be useful. 
The ratio 7(537.0 A)/J(584.3 A) would be expected to in- 
crease towards the edges of the network, but the variation 
would be relatively small, since the ratio of the enhancement 
factors of the two lines changes by only 25 per cent between 
log Te = 4.8 and log Te = 5.7. This factor is consistent with 
values of the ratio observed by MJ99, who found only small 
variations of the ratio, but a consistent increase in the ratio 
in regions of smaller absolute intensity. Enhanced excitation 
by non-local suprathermal electrons in the low TR would 
tend to produce a larger variation in the opposite sense with 
absolute intensities (Paper II). 

In a more refined approach to turbulent transport calcu- 
lations, improvements could be made to the simple geometry 
adopted here. Flux tubes with various inclinations could be 
modelled, allowing a more quantitative examination of how 
the spatial variation of helium emission may be affected. 
MJ99 found steeper temperature gradients in cell interiors 
than in the network, but these gradients may be across mag- 
netic structures. In the expanding network picture, the tem- 
perature gradient along the field is smaller, which would re- 
duce the enhancement factor, as would the higher pressures 
found by MJ99 in cell interiors. The greater helium line en- 
hancements required by MJ99's observations in cell interi- 
ors could possibly be caused by photon scattering from the 
boundaries; radiative transfer using two-component model 
atmospheres is needed to test this point. 

As described in Section [Tj], even significantly non- 
Maxwellian EVDFs are expected to have relatively little ef- 
fect on turbulent transport of helium in the quiet Sun. As 
enhancements of the helium line intensities by non-thermal 
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motions occur at temperatures much higher than the nor- 
mal temperature of peak line formation, and enhancement 
by non-local suprathermal electrons occurs largely at lower 
temperatures (Paper II), the two processes could be operat- 
ing simultaneously. If non-thermal motions do produce the 
enhancement factors calculated here, these would dominate 
over the effects of non-local electrons, provided departures 
from Maxwellian EVDFs are realistic (Paper II). 

The present calculations do not include radiative trans- 
fer, which could have significant effects on He I, but which 
would be less important for Hell. Radiative transfer calcu- 
lations including non-thermal motions (and possibly non- 
Maxwellian EVDFs as well) are needed to place our conclu- 
sions on a firmer footing. 
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